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Abstract 

Using data collected with the CLEO III detector at the Cornell Electron Storage Ring, we have 
compared proton, lambda (charged conjugate modes are implicit) and meson {(f) and /2(1270)) 
production observed in gluon fragmentation vs. quark fragmentation. Two studies have been 
conducted: in the first, we corroborrate previous per-event particle yields in T ggg vs. e"'"e~ — >■ 
qq. In the second, we compare particle production in the photon-tagged process T(1S) gg'j 
with that in e+e" — > qq'y events. For each particle, we determine the 'enhancement' ratio, defined 
as the ratio of particle yields per gluon fragmentation event compared to quark fragmentation 
event. Thus defined, an enhancement of 1.0 implies equal per-event production in both gluon and 
quark fragmentation. In the photon-tagged analysis, we find an enhancement of order unity for 
protons and approximately 1.5 for A's. This measured proton enhancement rate is supported by a 
study of baryon production in Xb,2 p + X relative to X6,i p + X. The production of mesons 
having masses of order 1 GeV and /2(1270)) are found to be also approximately the same in 
gg^ vs. qq^y fragmentation. Overall, per-event baryon production in two-gluon fragmentation is 
considerably smaller than that observed in three-gluon decays of the T(1S). Our results for baryon 
production arc inconsistent with the predictions of the JETSET (7.4) fragmentation model. The 
results presented in this document are preliminary. 



*Submitted to the 33'''^ International Conference on High Energy Physics, July 26 - August 2, 2006, Moscow 
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I. INTRODUCTION 



Understanding hadronization, or the process by which elementary partons (gluons and 
quarks) evolve into mesons and baryons, is complicated by its intrinsically non-perturbative 
nature. Due to the fact that gluons carry two color indices whereas quarks carry only 
one, the intrinsic gluon-gluon coupling strength (C^=3) is larger than the intrinsic quark- 
gluon coupling strength (Cir=4/3). Radiation of secondary and tertiary gluons is therefore 
expected to be more likely when hadronization is initiated by a gluon rather than by a 
quark. This results in a greater number of final state hadrons as well as a larger average 
transverse momentum in the former case compared to the latter case. In the limit — cxd, 
the ratio of the number of hadrons produced in gluon-initiated jets to the number of hadrons 
produced in quark-initiated jets is expected, in lowest order, to approach the well-known 
color-counting ratio C^/Ci?=9/4|.l]. 

Many experiments have searched for, and found, multiplicity and jet shape differences 
between quark and gluon fragmentation. At energies, qqg events are distinguished by 
their three-jet topology. Within such events, quark and gluon jets can be separated by a 
variety of techniques including vertex tagging. Because gluons rarely fragment into heavy 
quarks, they will produce jets that form a vertex at the e^e~ interaction point. Quark jets, 
to the contrary, tend to form a detached vertex when the jet contains a long-lived bottom 
or charm quark. For light-quark events with gluon radiation, however, the assignment of 
final state hadrons to the initial state partons is generally more ambiguous and often relies 
on Monte Carlo simulations to determine the fraction of times that an observed hadron is 
correctly traced to a primary parton. 

Within the limits of their precision, previous studies at SLD found inclusive production of 
pions, kaons and protons to be equivalent for gluon-tagged and quark-tagged jets|2]- OPAL 
has measured charm production at the level of (3.20 ± 0.21 ± 0.38)% in gluon jets[3, 0|, 
more than an order of magnitude smaller than the rate observed in quark jets at the . 
ALEPhQ and DELPHlj^ both measured inclusive bottom production in gluon-tagged jets 
to be 2 — 3 X 10"'^, again considerably smaller than that expected from charge counting 
in quark fragmentation. Most directly comparable to our current work, OPAL has also 
compared inclusive and A production in gluon- vs. quark-tagged jets, finding inclusive 
production ratios again consistent with unity (0.94 ± 0.07 ± 0.07 and L18 ± 0.01 ± 0.17, 
respectively) 13] . 

The 10 GeV center of mass energy range offers a unique opportunity to probe quark 
and gluon fragmentation effects, without relying on Monte Carlo simulation to associate 
the final state hadrons with an initial state parton. The decay T(15') gg^ allows one to 
compare the gg system in a gg^ event with the qq system in e'^e^ q(pi events. In both 
cases, the system recoiling against the photon consists (to lowest order) of hadrons that have 
evolved from either a two-gluon or a quark-antiquark system. The properties of the recoil 
systems can then be compared directly.^ Additionally, the radiative transitions from the 
radially excited T states to the orbitally excited Xb triplet offers an opportunity to further 



Although there may be gluon radiation from the initial partons, we do not distinguish such radiation 
explicitly in this analysis. Thus, the states that we are comparing are, strictly speaking, gg^ and qcpj 
to lowest-order only; additional gluon radiation, to which we are not experimentally sensitive, may be 
present in many of the events in our sample. Without the ability to adequately identify additional gluons, 
such higher-order radiative effects are therefore implicitly absorbed into the experimental measurement. 
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probe fragmentation differences between decays of the J=0 and J=2 Xb states (decaying 
predominantly to two gluons) vs. decays of the J=l state (decaying primarily to qq(g) final 
states, with the gluon nearly on-shell). Statistical correlations between transition photons 
with inclusive production of particular final-state particles (X) allows a measurement of the 
relative yields of gg X:qq{g)^ X. 

In this analysis, we focus on the relative production rates of baryons (protons and A's) 
and heavy mesons (0 and /2(1270)) in gluon vs. quark fragmentation (charge conjugation 
is implied). A previous study noted apparent enhancements in the production of (p, A and 
protons in three-gluon decays of the T(lS)j^, albeit with low statistics and at no more 
than 2-3(7 significance. That initial study also found approximately one unit larger charged 
multiplicity for three-gluon fragmentation of the T(1S) compared to qq fragmentation at a 
comparable center-of-mass energy. With the limited statistics at that time, the additional 
unit of multiplicity could entirely be accounted for by enhanced three-gluonic production of 
baryons. We now have sufficient statistics to re-measure the three-gluon particle production 
rates, and also to compare, for the first time, inclusive production in two-gluon fragmentation 
vs. inclusive production in three-gluon fragmentation. 

Since then, other experimental data on quark/gluon fragmentation differences in the 
^/s ~10 GeV energy regime have become available, including (per event): 



1. A suppression of production of open charm in gluonic decays of the T resonances: 
(T(1S)^ ggg-^D + X)/(e+e- ^ D + X) ~0.02[9]. 

2. The observation that fragmentation of the J=l state of the Xb triplet (decaying to 
qqg) results in charm production comparable to the underlying continuum; no such 
charm production is observed in the two-gluon decays of the J=0 or J=2 states [lol|. 

3. An enhancement in production of hidden charm in gluonic decays of the T resonances: 



(T(1S)^ ggg -> J/^ + X)/(e+e- ^ J/^ + X) >5.3 



11 



4. Production of deuterons from resonant 3-gluon decays of both the T(1S) and T(2S) 
at the level of 10~^; no corresponding production of deuterons is observed from the 
continuum|l3]. Per event enhancements are >10. 

5. Approximately equal production of rj' in gluonic decays of the T resonance compared 
to T decays via qS: (T(1S)^ ggg ^ rj' + X)/(T ^ qq ^ rj' + X) ~2/3, integrated 
over momentum ll 311. 



II. DETECTOR AND DATA SAMPLE 

The CLEO III detector [H O is a general purpose solenoidal magnet spectrometer and 
calorimeter. This system is very efficient (e >98%) for detecting tracks that have transverse 
momenta (px) relative to the beam axis greater than 200 MeV/c, and that are contained 
within the good fiducial volume of the drift chamber (|cos^^| <0.94, with 6 defined as the 
polar angle relative to the beam axis). Below this threshold, the charged particle detection 
efficiency in the fiducial volume decreases to approximately 90% at pt ~100 MeV/c. For 
Pt <100 MeV/c, the efficiency decreases roughly linearly to zero at a threshold of pt ~30 
MeV/c. Just within the solenoidal magnet coil is the electromagnetic calorimeter, consisting 
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of 7800 thallium doped Csl crystals. The central "barrel" region of the calorimeter covers 
about 75% of the solid angle and has an energy resolution of 



E 



(%) 



0.6 



+ 1.14-0.01E, 



(1) 



^0.73 



where E is the shower energy in GeV. This parameterization translates to an energy res- 
olution of about 2% at 2 GeV and 1.2% at 5 GeV. Two end-cap regions of the crystal 
calorimeter extend solid angle coverage to about 95% of 47r, although energy resolution is 
not as good as that of the barrel region. The tracking system, RICH particle identification 
system and calorimeter are all contained within the 1.5 Tesla superconducting coil. Flux 
return and tracking chambers used for muon detection are located immediately outside the 
coil and in the two end-cap regions. 

We use the CLEO-III data collected at the narrow T resonances as a source of 5^(77 events 
and data taken just below the resonances, as well as the below-4S continuum (y^=10.55 
GeV) as a source of 007 events. The 7 in our qq'-j sample results primarily from initial 
state radiation (ISR)[l^. We compare events for which the fractional photon energy = 
Ey/Ebeam are the same.^ 

A. Data Samples Used and Event Selection 

We use data taken at the T(1S), T(2S), T(3S), and T(4S) resonances. 

We impose event selection requirements identical to those used in our previous study of 
inclusive direct photon production in T decays^3]- Those cuts are designed primarily to 
supress backgrounds such as two-photon collisions, QED events (including tau pairs), and 
beam-gas and beam-wall collisions. 

Luminosity, event count, and photon yields (x^ >0.5) are given in Table HI 

B. Event Backgrounds 

To determine the characteristics of resonant gg^ events, we must subtract the background 
arising from non-resonant qq'y and e~^e~ TT'y events produced in continuum e'^e~ anni- 
hilations at y/s = Mx(is)- This is done by direct scaling of the event samples collected off- 
resonance on the nearby continuum. 

In order to isolate continuum qqj events, tt'j contamination must be explicitly sub- 
tracted, using a Monte Carlo simulation of tau pair events. We find that ttj events comprise 
about 5% of the qq'j data sample passing the event selection cuts[l3|- Beam gas and two 
photon backgrounds were investigated and found to be negligibly small. The photon-tagged 
sample can also be contaminated by cases where the high-energy photon is actually a 7r° 
daughter. Figure ^ illustrates the fraction of photons for the on-lS resonance and below- 
4S continuum, respectively, that are produced from neutral decays (including not only 7r° 
decays but also rj, rj\ and uj decays) as determined from Monte Carlo simulations. 



^ This deviates from the previous convention, for which the scahng variable was the invariant mass of the gg 
and qq systems recoiUng against the hard photon {MrecoiU defined by Mrecoii = \/4^^bcam(l ~ / Eheam)) ■ 
We take the difference between resuhs using these two conventions as a systematic error. 
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TABLE I: Summary of data used in analysis. For each data set, we track the number of photons 
per unit luminosity, as well as the total number of observed hadronic events per unit luminosity. 

HadEvts denotes the total number of events in each sample identified as hadronic by our event 
selection requirements. The number of photons having scaled momentum (relative to Eh^am) greater 
than 0.5 is presented in the last column. 
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FIG. 1: Monte Carlo prediction for fraction of photons produced by neutral decays (such as 7r°, 
ry, ry', and a;) on the IS resonance (left) and on the continuum below the 4S resonance (right). 
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Integrated over all tag photon momenta, 7r° contamination comprises a ~15% background 
to the direct-photon sample. Monte Carlo simulations also indicate that the vr^ contami- 
nation tends to cancel when we take ratios of resonant photon production to continuum 
photon production. 

1. Particle Identification 

Our photon and particle identification procedures are identical to those developed in . 
Photon candidates are selected from showers with widths and patterns of energy deposition 
consistent with that of a photon, as opposed to neutral hadrons (e.g. merged vr^'s, K^, 
neutrons, etc.). To ensure that the events are well-contained within the CLEO detector, 
we require | cos6'^| < 0.707 {6^ defined as before as the polar angle between the beam axis 
and the direct photon). For protons (and antiprotons), we require that charged tracks 
have specific ionizatation and also RICH information consistent with that expected for true 
protons. For momenta less than 1 GeV/c, we require that the track particle identification 
information be inconsistent (at the level of 2a) with that expected for true pions. We also 
require that the proton momentum be greater than 400 MeV, so as to suppress beam-wall 
and fake backgrounds. For reconstruction of (f) (and /2(1270)) from kaons (and pions) we 
require that pairs of opposite charge have particle identification information consistent with 
their assumed identities. For the case of the /2, we require the pion energy to be greater than 
500 MeV, to enhance the signal to noise ratio. Lambdas are identified using the standard 
CLEO algorithms for reconstruction of detached vertices. 

2. Backgrounds to the Proton Sample 

We use Monte Carlo simulations to assess fake proton backgrounds. Figure |21 illustrates 
proton fakes for a sample of below- IS Monte Carlo continuum simulations. The solid black 
curve shows the number of all particles identified as protons that were also tagged as true 
protons. The red dashed (blue dotted, magenta dash-dot) curve corresponds to those par- 
ticles that were identified as protons, but that were tagged at the generator level as true 
kaons (pions, positrons). Proton backgrounds are observed to be present at the ~10% level 
and are expected to largely cancel in the ratio of ratios. 

Note that, for all proton and antiproton identification, we require that the minimum 
(anti-)proton momentum exceed 400 MeV to eliminate concerns regarding protons ranging 
out in the beampipe. 

C. Signal Definition 

In this analysis we measure particle enhancements in both the ggg and gg7 decays of the 
Upsilon system, relative to qq(7) production on the underlying continuum. Our definition of 
enhancement is given quantitatively as the continuum-subtracted resonance yield relative to 
the continuum yield. Thus defined, an enhancement of 1 indicates that a given particle is pro- 
duced as often (per event) on the continuum as on the resonance. Note that our definition of 
'continuum' here means both below-resonance continuum, as well as resonance— ^qq through 
vacuum-polarization; i.e., all e^e~ — >qq-like processes which must be explicitly subtracted 
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FIG. 2: (Left) Proton fakes for a sample of below-lS Monte Carlo simulations. The solid black 
curve shows the number of all particles identified as protons that were also tagged as true protons. 
The red dashed (blue dotted, magenta dash-dot) curve corresponds to those particles that were 
identified as protons, but that were actually kaons (pions, positrons). (Right) Same for on-lS event 
simulations. 

in determining the characteristics of 3-gluon resonant decays.^ Furthermore, note that for 
the T(2S) and T(3S) data, there is no subtraction of cascades to lower T states or xt decays. 
Therefore, in what follows, "T(2S)" denotes a sum over T(2S) direct, T{2S) — > T(1S)+X 
and T(2S)— > xt- Assuming the direct decays of the T resonances are identical, an T(2S) 
enhancement smaller than that of the T(1S) implies that the enhancements from the first 
and third processes enumerated above are therefore smaller than for the T(1S). 

In general we have two continuum-subtraction options: we may determine enhancements 
for all resonances relative to the below-4S continuum (for which the statistics are largest, 
but the extrapolation in energy and run period is also largest) or we may find enhancements 
relative to their individual below-resonant continua using data not more than 50 MeV away 
from the resonance of interest. For mass- fitted particles we normalize exclusively to the 
below-4S contiuum, as the individual continuua (below-lS, -2S, and -3S) do not have high 
enough statistics to yield well-fitted mass peaks. For particle counts determined by the 
momentum spectra (protons and antiprotons), we normalize to both the below-4S continuum 
as well as the resonance-specific continuua and present the differences in the enhancements 
produced by the two as a systematic error. 



^ Vacuum polarization processes are subtracted by direct scaling of the continuum using the T — > 7* — > 
values tabulated previously 0|. 
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D. Particle production in three-gluon vs. qq events 

The previous CLEO-I analysis already observed significant enhancements of protons and 
lambda's produced in 3-gluon decays of the IS relative to the below-4S continuum |8|]. We 
repeat that analysis with our enlarged data set, as described below. 

Errors on particle yields are obtained from the error on the fit if the particle count is 
obtained by fitting a mass peak (A, 0, or by the square root of the total count if the 
particle count is obtained from the momentum spectrum (p,p). For the ggg analysis, we 
determine enhancements as a function of scaled momentum and also calculate momentum- 
integrated enhancements for each particle, to allow comparison with previous results. 

E. gg7 Analysis 

For the gg7 analysis we normalize the particle counts in a given bin to the photon counts 
in the same bin. In this case the bins are determined by the photon momentum and not 
the particle momentum (as is the case in the ggg analysis). For each bin, we then find the 
fractional contamination of resonance photons due to the underlying continuum. Once this 
fractional contamination is known, the resonance yield can be extracted by straight-forward 
algebra. As a check of the numerical procedure, we have verified that the enhancements 
obtained for T(4S) (which should produce no photons with >0.5, since its width is 
dominated by BB) relative to the below-4S continuum are consistent with zero. 



III. RESULTS 

A. ggg Enhancements 

1. Baryon Enhancements 

Figure 01 presents our A enhancements binned according to scaled momentum, that is the 
momentum of the particle divided by the beam energy. In the figure, blue square (green 
triangle, red diamond) symbols correspond to enhancements on the IS (2S, 3S) resonance. 
Closed symbols are data and open symbols are JETSET 7.4|3| event generator simulations 
followed by the full CLEOIII GEANT-based[ii] Monte Carlo detector simulation. We note 
that the acceptance in the lowest momentum bin is very sensitive to absorption of protons 
in beam-pipe. 

Figure |3] shows the proton and antiproton enhancements. The consistency between the 
two indicates that beam-wall and beam-gas backgrounds are not substantial. 

2. 4> and f 2 (1270) Enhancements 

Figure El shows (j) enhancement results binned according to scaled momentum. Symobls 
are as above with blue square (green triangle, red diamond) symbols correspond to enhance- 




.Continuum 



2 




9 



ggg A Enhancement 



■ 1S Data Corrected 
A 2S Data Corrected 
▼ 3S Data Corrected 
□ 1S Generator MC 

■ 2S Generator MC 
V 3S Generator MC 



Y 

m 



0.00 



4 



-4 — 4- 



0.25 0.50 0.75 

Scaled A Momentum 



1.00 



FIG. 3: Raw enhancements for ggg— > A + X binned according to scaled momentum {pa/ Ebeam)- 
Blue square (green triangle, red diamond) symbols correspond to enhancements on the IS (2S, 3S) 
resonance. Closed symbols are data, open symbols are JETSET Monte Carlo. 
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FIG. 4: (Left) Raw enhancements for ggg p + X binned according to scaled momentum. Statis- 
tics of continuum subtraction are largest in lower momentum bins. Blue square (green triangle, 
red diamond) symbols correspond to enhancements on the IS (2S, 3S) resonance. Closed symbols 
are data, open symbols are JETSET Monte Carlo. (Right) Same for antiprotons. 
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FIG. 5: (Left) Raw enhancements for ggg (j) binned according to scaled momentum. Blue square 
(green triangle, red diamond) symbols correspond to enhancements on the IS (2S, 3S) resonance. 
Closed symbols are data, open symbols are JETSET Monte Carlo. (Right) Enhancements for 
/2(1270). 

merits on the IS (2S, 3S) resonance. Closed symbols are data and open symbols are JETSET 
Monte Carlo. Here, we have normalized the resonant production at 9.46 GeV to the contin- 
uum production at 10.55 GeV. Therefore, for production, the lowest momentum bins for 
the resonance are particularly sensitive to low-momentum kaon acceptance. Figure |S1 also 
shows the f2 enhancement results binned according to scaled momentum. The f2 peak is not 
well-defined at low momentum (lowest two bins). No Monte Carlo comparison is presented 
since our current Monte Carlo, by default, will not generate tensor particles. 



3. Particle Momentum- Integrated Enhancements 

Figure IHl shows the particle momentum-integrated enhancements for each particle. We 
note that the baryons (A, p, p) have enhancements greater than 1, while the mesons (0, 12) 
have enhancements less than 1. Our results are, in general, numerically consistent with the 
prior CLEO-I analysis, albeit with considerably higher statistics. 



B. gg7 Enhancements 

There are sufficient statistics to present binned enhancements (in this case binning is done 
according to tagged photon momentum, not according to individual particle momentum) for 
A, protons and p. For all particles, we also present momentum-integrated enhancements. 
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FIG. 6: Compilation of momentum-integrated enhancements for ggg events. Blue square (green 
triangle, red diamond) symbols correspond to enhancements on the IS (2S, 3S) resonance. Closed 
symbols are data, open symbols are JETSET Monte Carlo. The x-axis is arbitrary, chosen to 
display all particles on one plot: A is plotted at x=l, protons at x=2, p at x=3, (j) x=4, and 
f2 at x=5. Systematic errors and relative efficiencies have now been included for this compilation. 
The CLE084 study did not measure an enhancement for /2(1270) and also only presented a single 
enhancement for the sum of protons and antiprotons. 



1. Baryon Enhancements 

Figure [7| shows A results binned according to scaled photon momentum, with the color 
scheme as before. Errors are large in the highest bin due to limited statistics. 

Figure |H1 shows proton and antiproton enhancement results binned according to scaled 
photon momentum. We immediately note somewhat smaller enhancements compared to the 
3-gluon enhancements presented previously. 
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FIG. 7: Raw enhancements for ggj A + X binned according to scaled photon momentum. 
We have integrated over all A momenta to make these plots, in events having a photon with the 
indicated x-value. Blue square (green triangle, red diamond) symbols correspond to enhancements 
on the IS (2S, 3S) resonance. Closed symbols are data, open symbols are JETSET Monte Carlo. 

2. Momentum- Integrated Enhancements 

Figure El shows the momentum- integrated enhancements for each particle. We note that 
the enhancements are considerably mitigated in comparison with the case of three-gluon 
fragmentation. 



IV. Xb,j^P+^ 

Photon transitions of the T(2S) and T(3S) to the Xb states allow us to measure the baryon 
yields, in association with a radiative transition photon 'tag'. Typical photon tag energies in 
this case are of order 80-160 MeV. Due to the large background in this region and therefore 
limited statistical power, the data permit only an extraction of the proton and antiproton 
enhancements. Of particular interest is the proton yield in Xb.2 vs. Xb,i decays; the former is 
expected to be dominated by decays via two-gluons, the latter is expected to be dominated 
by decays to qq(g), with the gluon expected to carry away very little momentum. 

We first conduct a Monte Carlo study to determine the relative efficiency of reconstructing 
a J=2 transition photon relative to J=l event, and also the efficiency when we require 
that a proton be found in addition to the transition photon. We compile statistics on the 
Xb{') —>■ p{p) + X analyses, separately for J=0/J=1 and for J=2/J=l. For the latter, we 
fit the signal to a double Gaussian, so that the errors in the areas of the two peaks are 
correlated and the statistical significance of the signals enhanced. For the former, we simply 
fit two Gaussians directly. We find that the efficiency for reconstructing a photon-proton 
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FIG. 8: (Left) Raw enhancements for gg7 p binned according to photon momentum. We have 
integrated over all proton momenta to make these plots, in events having a photon with the 
indicated x-value. Blue square (green triangle, red diamond) symbols correspond to enhancements 
on the IS (2S, 3S) resonance. Closed symbols are data, open symbols are JETSET Monte Carlo. 
(Right) Same for antiprotons. 

correlation in a two-gluon decay is approximately 95% that for a photon-proton correlation 
in qq(g) event. 

To check the sensitivity to our particle identification criteria, we have compared results 
using very tight proton identification requirements (with a reduction in efficiency by more 
than 50%) vs. the 'standard' loose proton identification criteria used above. We obtain a 
comparable correction factor using more restrictive particle identification criteria. 

Results are presented in the Table ITTl We note that the observed enhancements are, again, 
smaller than those observed in comparing three-gluon fragmentation from the T resonance 
with qq fragmentation. 



V. CROSS-CHECKS AND SYSTEMATICS 

In order to verify our procedures and probe possible systematic uncertainties, two pri- 
mary cross-checks were employed. We first compare the Monte Carlo enhancements at the 
event generator-level with those determined after the generated events are processed through 
the full CLEO-III detector simulation ("detector-level"), as a function of momentum. In 
general, these enhancements will vary for several reasons, including differences in: a) the 
efficiencies for finding recoil particles in qq7 vs. gg7 events resulting from angular distribu- 
tion, event multiplicity, and particle momentum differences, b) event selection efficiencies, 
c) vr" contamination levels, and d) recoil center-of-mass discrepancies, e.g., between the con- 
tinuum data under the T(1S) resonance vs. the below-T(4S) continuum. In cases where the 
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FIG. 9: Momentum-integrated enhancements for gg7. Blue square (green triangle, red diamond) 
symbols correspond to enhancements on the IS (2S, 3S) resonance. Closed symbols are data, open 
symbols are JETSET Monte Carlo. The x-axis is arbitrary, chosen to display all particles on one 
plot: A is plotted at x=l, p at x=2, p at x=3, (j) x=4, and £2 at x=5. Systematic errors have 
been included and relative efficiency corrections have been applied. 

generator-level and detector-level enhancements are statistically inconsistent with each at 
the 2a level, we use the ratio between the generator- level and detector- level enhancements 
as a correction factor and conservatively take half of the amount by which this correction 
deviates from unity as an estimated systematic error. (Note that these corrections have 
already been incorporated into the results presented in Figures and IHl and E)) . Figures ITUl 
and ^2 shows the comparison of proton enhancements determined at the event-generator 
vs. post-detector-simulation levels of Monte Carlo simulation. Integrated over momentum, 
typical corrections are typically of order 10%. 

To test the sensitivity of our analysis procedures across different running periods, 
we have calculated the enhancements for photon-tagged T(4S) on-resonance events vs. 
photon-tagged below-T(4S) continuum events, spanning the full CLEO-III data set. Since 
T(4S)— s>BB~100%, we expect that any event having a photon with > 0.5 is a continuum 
event. Hence, the calculated enhancement should be zero. In all cases, save for antipro- 
tons, we find good agreeement between the below-4S continuum particle yields per photon 
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FIG. 10: (Left) Scaled momentum binned enhancements for ggg p + X. Blue square (green 
triangle, red diamond) symbols correspond to enhancements on the IS (2S, 3S) resonance. Closed 
(open) symbols are generator (detector) level Monte Carlo enhancements. (Right) Same for an- 
tiprotons. 
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FIG. 11: (Left) Scaled momentum binned enhancements for gg7 decays to p. Blue square (green 
triangle, red diamond) symbols correspond to enhancements on the IS (2S, 3S) resonance. Closed 
(open) symbols are generator (detector) level Monte Carlo enhancements. (Right) Same for an- 
tiprotons. 
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TABLE II: Summary of inclusive proton (and antiproton) results for Xb,j decays. For checks 
of internal consistency, data have been separated into sub-samples, labeled with capital Roman 
letters. For J=2 relative to J=l, e.g., the scale of systematic uncertainties is set by the constancy 
of the value across datasets (r.m.s.~0.03), the magnitude of relative efficiency corrections (~0.05) 
and the consistency of results obtained using different particle identification criteria. 
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tag, and the on-4S particle yields per photon tag. For antiprotons, we find deviations from 
the null expectation at the level of ~10%, and incorporate these deviations (bin-by-bin in 
momentum) into our total systematic error for that particular case. 

VI. DISCUSSION AND SUMMARY 

We compare, for the first time, particle production in two-gluon vs. quark- ant iquark 
fragmentation. We find that, in particular, baryon production (per event) in two-gluon 
decays is somewhat smaller than that observed in three-gluon decays, particularly in the 
case of A production. This result is qualitatively at variance with the conclusion that, 
e.g., the thrust and charged multiplicity distributions of Xb,o and Xh,2 two-gluon decays 
agreed well with T(1S)— ggg and that the thrust and charged multiplicity distribution of 
qqg agreed with continuum e"'"e~ — > qq events jiO]. This effect is not reproduced in the 
current JETSET 7.4 Monte Carlo event generator. For protons, which represent our highest- 
statistics sample, our results are inconsistent with a model where baryon production in gluon 
fragmentation is only a function of the available center-of-mass energy. Further detailed 
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comparisons with models based on either string or independent parton fragmentation may 
help clarify the production mechanisms. 
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